Paths to Synchronization on Complex Networks 
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The understanding of emergent collective phenomena in natural and social systems has driven the interest 
of scientists from different disciplines during decades. Among these phenomena, the synchronization of a set 
of interacting individuals or units has been intensively studied because of its ubiquity in the natural world. In 
this paper, we show how for fixed coupling strengths local patterns of synchronization emerge differently in ho- 
mogeneous and heterogeneous complex networks, driving the process towards a certain global synchronization 
degree following different paths. The dependence of the dynamics on the coupling strength and on the topology 
is unveiled. This study provides a new perspective and tools to understand this emerging phenomena. 
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In 1998 Watts and Strogatz in an effort to understand the 
synchronization of cricket chirps, which show a high degree of 
coordination over long distances as though the insects where 
"invisibly" connected, end up with a seminal paper about the 
small-world effect (H that was the seed of the modern theory 
of complex networks El 01 ■ Many natural and man-made 
networks have been, since then, successfully described within 
this framework. Nevertheless, the understanding of the syn- 
chronization dynamics in complex networks remains a chal- 
lenge. 

The synchronization of non-identical interacting units oc- 
cupies a privileged position among emergent collective phe- 
nomena because of its various applications in Neuroscience, 
Ecology, Earth Science, among others JH S 0]- One of the 
most successful attempts to understand it is due to Kuramoto 
IH 01, who analyzed a model of phase oscillators coupled 
through the sine of their phase differences. The Kuramoto 
model (KM) consists of a population of N coupled phase os- 
cillators where the phase of the i-th unit, denoted by 6i(t), 
evolves in time according to 
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where stands for its natural frequency, Ay is the cou- 
pling strength between units and Ay is the connectivity ma- 
trix (Ay = 1 if i is linked to j and otherwise). The original 
model studied by Kuramoto assumed mean-field interactions 
with Aij = l,Vi ^ 3 (all-to-all) and Ay = JC/N, Mi, j. The 
model can be solved in terms of an order parameter r that 
measures the extent of synchronization in a system of N os- 
cillators as: 
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where ^ represents an average phase of the system. The pa- 
rameter < r < 1 displays a second order phase transition in 
the coupling strength, being r = the value of the incoherent 
solution, and r = 1 the value for total synchronization. 



The synchronization problem has been solved in some other 
cases, mainly those where a mean-field approach is also valid 
lioll . Unfortunately, the mean-field approach requires of sev- 
eral constraints that are not usually fulfilled in real systems. 
Natural, social and technological systems show intricate pat- 
terns of connectivity between their units that are, nowadays, 
described as complex networks JHQ]. The problem of syn- 
chronization in complex networks inherits the technical diffi- 
culties of the non mean-field approaches and incorporates new 
questions to be considered: What are the new pertinent param- 
eters to deal with synchronization? and, What is the role of the 
topology in the synchronization process? Several works have 
partially addressed these issues by stu dying th e stability of the 
synchronized state 1 11 , 12 , T^l using the Mas- 

ter Stability Function (MSF) formalism 111 81 . However, the 
onset of synchronization, which posses more theoretical and 
phenomenological challenges is much less explored, and only 
a few works have dealt with the study of the whole synchro- 
nization dynamics in specific scenarios lfl9l 2oL 21. 22, 231. 

The main goal of this Letter is to study the synchroniz- 
ability of complex networks as a function of the coupling 
strengths. To do this, first, we propose and discuss a new mea- 
sure of synchronization for the KM in complex networks. Sec- 
ond, we scrutinize and compare the synchronization patterns 
in Erdos-Renyi (ER) and scale-free (SF) networks and show 
that even in the incoherent solution, r = 0, the system self- 
organizes towards synchronization following different paths. 
Our study reveals that the synchronizability of these networks 
does depend on the coupling between units, and hence, that 
general statements about their synchronizability are eventu- 
ally misleading. Moreover, we show that even in the inco- 
herent solution, r = 0, the system is self-organizing towards 
synchronization. 

Let us start by mapping the KM model in finite complex 
networks as 

d6 N 
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where A is a constant. 
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FIG. 1: (color online) Evolution of a, the KM order parameter de- 
fined in Eq. (|2), and b the fraction of synchronized links ru n ^, Eq. 
j4j, as a function of A. The curves separate when the incoherent so- 
lution for SF networks destabilizes. The figure clearly illustrates that 
the synchronizability of the networks does depend on the value of the 
coupling strength. Both plots are represented for Erdos-Renyi (ER) 
and scale-free (SF) networks as indicated. The size of the networks 
is TV = 1, 000 and their average degree is (k) — 6. The exponent of 
the SF network is 7 = —3. 



We study the dynamics of Eq.© in ER and SF networks, 
preserving the total number of links, Ni and nodes, N for 
a proper comparison 112411 . We concentrate in two aspects: 
global and local synchronization. First, we follow the evo- 
lution of the order parameter r, as A increases, to capture the 
global coherence of the synchronization in the networks. Sec- 
ondly, we propose and follow the same evolution for a new pa- 
rameter, ru n k . This parameter measures the local construction 
of the synchronization patterns and allows for the exploration 
of how global synchronization is achieved. We define 
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that represents the fraction of all possible links that are syn- 
chronized in the network (averaged over a large enough time 
interval At, after the system relaxes at some large time t r ), 
being T; the set of neighbors of node i. 

We solved Eq.© using a 4 th order Runge-Kutta method 
for different values of A, with a uniform distribution of nat- 
ural frequencies g(cu) in the interval [— n, ir] up to achieving 
the stationary state. The networks are built following a model 
1 2511 that generates a one parameter family of complex net- 
works. This parameter, a G [0, 1], measures the degree of 
heterogeneity of the final networks. A network of size N is 
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FIG. 2: (color online) Size of largest synchronized connected com- 
ponent (GC) and number of synchronized connected components 
(Nc), as a function of A for the different topologies considered. De- 
spite r being vanishing and hence no global synchronization is yet 
attained, a significant number of clusters show up. This indicates 
that for any A > the system self-organizes towards macroscopic 
synchronization. The network parameters are as in Fig.Q] 



generated starting from a fully connected core of mo nodes 
and a set U(0) of N — mo unconnected nodes. At each time 
step, a new node (not selected before) is chosen from U(0) 
and linked to m other nodes. Each of the m edges is linked 
with probability a to a randomly chosen node (avoiding mul- 
tiple and self-connections) from the whole set of A — 1 re- 
maining nodes and with probability (1 — a) following a linear 
preferential attachment strategy 12611 . Repeating these steps 
(A — mo) times, networks interpolating between the limiting 
cases of ER (a = 1) and SF (a = 0) topologies are generated 
& 

In Fig. Q] we represent the evolution of both order parame- 
ters, r and ru n k, as a function of the coupling strength A. The 
global coherence of the synchronized state, represented by r, 
shows that the onset of synchronization first occurs for SF net- 
works. A detailed finite size scaling analysis performed for 
both topologies shows that the critical value of the effective 
coupling, A c , corresponds in SF networks to Af F = 0.05(1), 
and in ER networks to Xf R — 0.122(2), accordingly with 
Fig.Q] If A is further increased, there is a value at which r for 
the ER crosses over the SF curve. From this value up in A, the 
ER network remains slightly more synchronized than the SF 
network. 

The behavior of ru n k shows a change in synchronizability 
between ER and SF and provides additional information. In- 
terestingly, the nonzero values of ru n k for A < A c indicate 
the existence of some local synchronization patterns even in 
the regime of global incoherence (r rj 0). Right at the onset 
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FIG. 3: (color online) Synchronized clusters for several values of A 
for the two different topologies studied (ER and SF). These networks 
are made up of 100 nodes, in order to have a sizeable picture of the 
system. The evolution of local synchronization patterns is always 
agglomerative, however, it follows two different routes. In the ER 
case, the growth of the GC proceeds by aggregation of small clusters 
of synchronized nodes, while for the SF network the central core 
groups the smaller clusters around it. 



of synchronization for the SF network, its ru n k value deviates 
from that of the ER. While the synchronization patterns con- 
tinue to grow for the ER network at the same rate, the forma- 
tion of locally synchronized structures occurs at a faster rate in 
the SF network. Finally, when the incoherent solution in the 
ER network destabilizes, the growing in its synchronization 
pattern increases drastically up to values of ru n k comparable 
to those obtained in SF networks and even higher. 

The above results show that statements about synchroniz- 
ability are dependent on the coupling strength value. Addi- 
tionally, the previous discussion suggests that synchronization 
is attained following two different paths that depend on the un- 
derlying topology. We have studied the characteristics of the 
synchronization patterns along the evolution of ru n k . Syn- 
chronization patterns are formed by pairs of oscillators, physi- 
cally connected, whose phase difference in the stationary state 
tends to zero. Note that the contribution into Eq.© of every 
pair of connected oscillators can be written in terms of a ma- 
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matrix is filtered using a threshold T such that the fraction 
of synchronized pairs equals rimk ■ In this wav > if f ij > T 
oscillators i and j are considered synchronized and the syn- 
chronized patterns are extracted. 

In Fig. |2] we represent the number of synchronized clus- 
ters and the size of the largest one (GC) as a function of A. 
The local information extracted from it is unveiling an aston- 
ishing and novel feature of the synchronization process that 
can not be derived from Fig. [T] and that in some sense is 
counterintuitive. The emergence of clusters of synchronized 
pairs of oscillators (links) in the networks shows that for val- 
ues of A < Af F , i.e., still in the incoherent solution r = 0, 
both kind of networks have developed a largest cluster of syn- 
chronized pairs of oscillators involving 50% of the nodes of 
the network, and an equal number of smaller synchronization 
clusters. From this point on, in the SF network the GC grows 
and the number of smaller clusters goes down, whereas for 
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FIG. 4: (color online) The plot shows the correlation between the 
likelihood that a node belongs to the GC of pairs of synchronized 
oscillators and its degree k as a function of the coupling strength A. 
This probability ,P GC (k), is color-coded as indicated in the right 
panel. The figure convincingly demonstrates that highly connected 
nodes are those that recruit poorly connected nodes as the GC grows. 
Network parameters are those used in Fig.[T] 



the ER network the growth exploits. These results indicate 
that although SF networks present more coherence in terms 
of r and ru n k, the microscopic evolution of the synchroniza- 
tion patterns is faster in ER networks, being these networks 
far more locally synchronizable than the heterogeneous ones. 

The observed differences in the behavior at a local scale 
are rooted in the growth of the GC. It turns out that for the 
ER networks, many different clusters of synchronized pairs 
of oscillators merge together to form a GC when the effective 
coupling is increased. The coalescence of many small clus- 
ters leads to a giant component of synchronized pairs that is 
almost the size of the system once the incoherent state desta- 
bilizes. This is not anymore the case for SF networks, where 
oscillators are incorporated to the GC practically one-by-one 
(forming new pairs) in terms of A (or ru n k), but starting from 
a core made up of half the nodes of the network. This picture 
is confirmed in Fig. [3] where we have represented the evolu- 
tion of local synchronization patterns in ER and SF networks 
for several values of A. The ultimate reason behind these two 
different routes to complete synchronization is the heteroge- 
neous character of the SF network and the role played by the 
hubs. In Fig. |H we have plotted the probability that a node 
with degree k belongs to the GC as a function of its degree 
k and the coupling A for the SF network. This probability is 
an increasing function of k for every A, hence the more con- 
nected a node is, the more likely it takes part in the cluster 
of synchronized links. Recently 12311 . Zhou and Kurths have 
reported the study of hierarchical organization in complex net- 
works, using the MSF and a mean-field approach in the weak 
coupling limit. Our results thus substantiate and generalize 
those about the role of hubs in the synchronization process 
presented in ll23"ll . 
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In summary, we have shown that synchronizability of com- 
plex networks is dependent on the effective coupling A among 
oscillators. For small values of A, SF networks outperform 
ER topologies, but the tendency is reverted for intermediate to 
large values of the coupling. On the other hand, the detailed 
analysis of evolution of patterns of synchronization showed 
that there are two radically different mechanisms to attain syn- 
chronization. In the presence of hubs, a giant component of 
synchronized pairs of oscillators forms and grows by recruit- 
ing nodes linked to them. On the contrary, in homogeneous 
structures, many small clusters first appear and then group to- 
gether through a sharp merging process. These results are as 
far reaching as the ones obtained for percolation and epidemic 
spreading on top of homogeneous or heterogeneous graphs, 
where the radical differences of the system's dynamics are 
rooted in the topology of the underlying networks, demon- 
strating that the same behavior may hold for nonlinear dy- 
namical systems coupled to complex structures. However, at 
variance with percolation processes, here the synchronization 
patterns could be directly related to the growth and evolution 
of the network. Therefore, we have naturally incorporated a 
dynamics relevant for the emergence of cooperative behavior, 
showing that the same organizing principle may drive network 
evolution, i.e., if synchronization is a relevant issue, natural 
networks (whether they are homogeneous or heterogeneous 
in degree) can indeed be efficient by adaptively selecting their 
coupling strengths. Our study then opens new paths to clar- 
ify how synchronization is attained in complex topologies and 
provides new tools to analyze this ubiquitous phenomenon. 
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